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Abstract 

New integrable multi-atom matter-radiation models with and without rotating wave approx- 
imation (RWA) are constructed and exactly solved through algebraic Bethe ansatz. The models 
with RWA are generated through ancestor model approach in an unified way. The rational case 
yields the standard type of matter-radiaton models, while the trigonometric case corresponds to 
their q-deformations. The models without RWA are obtained from the elliptic case at the Gaudin 
and high spin limit. 

42.50 Pq, 03.65 Fd, 32.80 -t 

1 Introduction 

It is amazing to note that real systems, like those in quantum optics induced by resonance interaction 
between atom and a quantized laser field, in cavity QED both in microwave and optical domain 
[1, 2], in trapped ion interacting with its center of mass motion irradiated by a laser beam [3, 4] 
etc., all involving complex matter radiation (MR) interactions can be described so successfully by 
simple models, like Jaynes-Cummings (JC) [5] and Buck Sukumar (BS) [6] models. Many theoretical 
predictions based on these models, like vacuum Rabi splitting , Rabi oscillation and its quantum 
collapse and revival etc. have been verified in maser and laser experiments. 

In such simplest interacting MR models the quantized radiation field is taken in single bosonic 
mode: h = be~ % ^5 1 + b'e^f*, while the atom is considered to be a two-level spin system with polar- 
ization vector S = a~e~ tuJat + <j + e +lulat . The interaction (h ■ S) therefore contains in general both 
fast (with frequency Wf + w a ) and slow oscillating (with frequency wj — w a ) components. However 
it is customary to neglect the fast oscillating part by considering approximation \wf — w a \ « Wf, 
which corresponds to the rotating wave approximation (RWA) and make the model solvable. This 
yields with the addition of free field and atomic excitation terms the well known JC model 

Hj C = uj f b ] b + uj a a z + a{b ] a~ + ba + ). (1) 
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However this approximation, which is justified only near the resonance point Uf ~ cj a , should be 
avoided in the general case [11], when additional counter rotating wave (CRW) terms must appear: 

H crw = (3(tfa + + ba-) (2) 

Moreover, for describing physical situations more accurately, one has to look for generalizations of 
the basic models, i.e. has to consider their q-deformations [7, 8], introduce higher nonlinearities 
[4, 9], involve multiple atoms [2, 9, 10] etc. However, while the exact solutions for the JC and the BS 
models together with their simple multi-atom extensions are known [12, 13, 14], the same appears to 
be no longer true for most of the above generalizations. More precisely, integrable models without 
RWA as well as those with explicit inter-atomic couplings are not known for most of the MR models, 
except for a few recent attempts [15, 16]. Moreover, though q-deformation, signifying introduction of 
anisotropy together with specific nonlinearity into the system, was considered earlier for a few models 
[7, 8], their multi-atom and integrable variants are not known. Therefore, it is indeed a challenge to 
find a scheme for generating integrable MR models having the desired properties mentioned above. 

2 Major aims and strategy of our construction 

To meet this challenge we construct new classes of integrable MR models with and without RWA. The 
result concerning our models with RWA has been reported recently [17], where we propose a general 
integrable system based on our ancestor Lax operator approach [18]. Through various reductions 
we generate from it a series of integrable MR models with explicit inter-atomic interactions. This 
includes such new multi-atom generalizations of JC, BS and trapped ion (TI) models as well as their 
integrable q-deformations involving quantum group. Moreover, since the construction of our models 
is based on a general Yang-Baxter (YB) algebra, we can also solve them in a unified way through 
the Bethe ansatz (BA). Our aim here is to review briefly this result and then present our new result 
obtained without RWA. 

Our strategy for constructing models with RWA is to start with a combination of Lax operators: 
T(A) = L S (X) T\f a Lj(A), where L S (X) involving bosons represents the radiation or the oscillatory 
mode and Lj(X) involving spin operators represents iV a -number of atoms. By construction it must 
satisfy the quantum YB equation (QYBE) R(X - n)T(X) <g> T(jj) = (I <g> T(jj))(T(X) <8> I)R(X - n), 
yielding the integrability condition [t(A),t(/z)] = 0, with mutually commuting set of conserved 
operators, obtained through expansion r(A) = trT(X) = J2a^a^ a [19]. For constructing the bosonic 
Lax operator we implement our ancestor model approach [18] and from a generalized Lax operator 
linked with a generalized quadratic algebra or its quantum deformation generate different L S (X) 
through various bosonic realizations. Our standard MR models are linked to the Lie algebra and the 
rational .R-matrix of the xxx spin chain [19], while their (/-deformations are related to the quantum 
algebra and the trigonometric i?-matrix of the xxz chain [20]. 

For constructing integrable MR models without RWA, however, as we explain below, the above 
approach fails due to the difficulty in bosonizing the Sklyanin algebra . Therefore for such models, 
which are linked to the elliptic i?-matrix of the xyz spin chain, we have to change our strategy and 
switch over to the Gaudin limit together with a high spin limit. 
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3 Integrable MR models with RWA 



We concentrate first on standard MR models and recall that in the rational case the 2x2 ancestor 
Lax operator may be given as [18] 



c?(A + s 3 )+cl, 8~ 
s+, c°(A- S 3 )-c 2 



with operators s satisfying a quadratic algebra 

[s + ,s~] = 2m + s 3 + m~, [s 3 ,s ± ]=±s ± , [m ± ,-]=0. (4) 

The central elements are expressed through arbitrary parameters appearing in (3) as m + = 
c°C2, m~ = c\c2 + c\c\ and as it is easy to see, their different choices reduce (4) to different simple 
algebras: 

i) su(u), at m + = 1, mT = 0, ii) su(l, 1), at m + = —l,m~ = 0, 
hi) bosonic, at m + = 0,m~ = —1, iv) canonical, at m + = m~ = (5) 

and the corresponding limits yield from (3) the respective Lax operators. In case i), choosing c„ = 1 
and c\ = — Crj = Cj, (3) we get the spin Lax operator Lj(A),j £ [l,iV ] describing N a atoms with 
inhomogeneity parameters Cj , while the other cases can reduce to different types of bosonic Lax 
operators linked to the radiation mode. 

Remarkably, such reductions yield in a unified way new integrable multi-atom BS, JC and TI 
models, at the limits ii), hi) and iv) of (5). Thus case ii) with choice c° = — cfj = 1, c\ = c\ = c, yields 

N a N a 

H bs = oj f s 3 + J2 (uajS* + a(s + S- + 8-Sfj) +aJ^(S-S+ - S+Sr), (6) 

j i<3 

which with a bosonic realization of su(l, 1): s + = \/TVb^,s~ = b\fN,s 3 = N + ^ and taking 
the spin-s operator as S = -^Y^f&ki represents a new integrable multi-atom BS model with inter- 
atomic interactions and nondegenerate atomic frequencies. At iV a = 1 the matter-matter interactions 
obviously vanish, recovering the known model [14]. The radiation frequency ujf and the atomic 
frequencies u> a j in our models are defined in general through the inhomogeneity parameters as 

w / = Yl W J> w i = a ( c i - c °) c i> w «i =UJ f~ w i + a ( c i + 4) (7) 
j 

Similarly, under reduction hi), choosing c\ = a, c\ = 0,c\ = c,c\ = —a~ l and taking direct 
bosonic realization s~ = b,s + = b^s 3 = b^b, we obtain from the same (3) an integrable multi-atom 
JC model with matter-matter coupling. 

From reduction iv), by fixing the parameter values as c° = — 1, c\ = c, = c\ = and considering 
a consistent realization through canonical variables: = e^ lx , s 3 = p + x, we generate an integrable 
multi-atom TI model with full exponential nonlinearity: a(e~ tx S + + e lx S~). Such models however 
are difficult to solve by standard BA, due to the absence of pseudovaccum like in the Toda chain. 
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For constructing integrable q-deformed MR models the strategy is the same; only one has to start 
now from the trigonometric type ancestor Lax operator involving generalized q-deformed operators 
[18]. This yields an integrable system with Hamiltonian in the form 

H q MR = H d + (s+S~ + s~ S q ) sin a, 

H d = -ic cos(aA) + csin(aX), X = (s 3 q - S q + u), (8) 

which represent a new class of MR models with S q belonging to the quantum group U q (su(2)) and 
s q to a more general quantum algebra [18]. 

Integrable q-deformed BS model is obtained from (8) by realizing s q through q-oscillator: s+ = 

^J[N] q b q , s q = b q ^J[N] q , s q = N + 2, and using quantum spin operator S q as its co-product [20] : 

= J2jQ '^- Jk< i' Tk o-fq^- Jl> 3 (Tl , S z = At s = 1, one gets an integrable version of an earlier 

model [8]. 

Similarly realizing s q = b q , s~ = b q , s 3 = N yield a new integrable q-deformation of the JC model, 
and realization through canonical operators results an integrable q-deformed TI model. 

We emphasize again that all integrable MR models we proposed, similar to their unified construc- 
tion, allow exact BA solutions also in a unified way. More details on the above models and their 
exact solutions can be found in [17]. 

4 Integrable MR models without RWA 

The basic idea for constructing integrable matter-radiation models with RWA, as we have explained 
above, is to consider first the xxx or the xxz spin chain with arbitrary spins and then replace one of 
the spins by its bosonic realization. The first case belongs to the rational class, where the spins satisfy 
the Lie type algebra, while the second case falls in the trigonometric class, where the higher spin 
operators satisfy the quantum algebra. Fortunately both these algebras allow bosonic realization, 
which yields the required radiation mode. 

For constructing more general integrable models with additional CRW terms (like (2) and its 
multiatom generalizations), one may expect to apply the same idea, but for the xyz spin chain. 
However in this anisotropic case, linked to the elliptic i?-matrix, the higher spin operators satisfy the 
Sklyanin algebra [21], for which unfortunately no bosonic realization is known. In search for a way 
out, one observes that the Sklyanin algebra fortunately reduces again to the usual sit (2) algebra at the 
elliptic Gaudin limit [22], allowing the required bosonization. However, a direct bosonic realization of 
a single spin operator in Gaudin models fails again , since it results either nonintegrable or unphysical 
models. The situation is saved finally by a further limit of s — > oo, for the bosonic mode, which yields 
the desired integrable multi-atom JC type model with nontrivial CRW terms. We should mention 
that similar high spin limit has been used recently for constructing integrable multi-atom JC model, 
but with RWA [16]. 

4.1 Elliptic Gaudin model 

Since the first step in our construction is the elliptic Gaudin model (EGM), we review briefly the 
related results from [22]. EGM is obtained at the a — > limit from the inhomogeneous xyz spin 
model and therefore all relevant objects like the Lax operator, i?-matrix, conserved quantities, as 
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well as important formulas in the Bethe ansatz method, like the eigenfunctions and the eigenvalues 
together with the Bethe equations etc. concerning this model can be derived directly from those of 
the xyz model [23] at the said limit. The Lax operator with inhomogeneity parametrs z n and the 
r-matrix of the EGM can be expressed through elliptic functions in the form [22] 

3 3 

L n (A) = ]T w a (X - z n )S a n a\ r(A) = £ w a (X)a a ® a\ n = 0, 1, . . . , N a 

a=l a=l 

where 101(A) = —(A), w 2 (X) = —(A), w 3 (X) = —(A). (9) 
sn sn sn 

Operators (9) satisfy a semiclassical YBE: [L n (X) (g) L n (p)] = [r(A — p), (L n (X) I + I <8> L n (/j,)] , 
obtained at a — ► limit from the QYBE. Mutually commuting conserved quantities -fT„ of the EGM 
are obtained from the expansion coefficient of the xyz transfer matrix 

r xyz (X,a) = I + a 2 T (2 HX) + 0(a 3 ) (10) 



as 

r^(X^z n ) = H n = £ ™ a (z„ - z m )S^ (11) 

with X)„ #n = 0. 

It is most crucial for our purpose to note that, higher spin operators S® appearing in the Gaudin 
model (9), are reduced from the generators of the Sklyanin algebra to simply those of 

[S'n'^'m] = 2p5 nm S n , [S n , S^] = i^rem'S'n) (12) 

with p = ± corresponding to su{2) (or su(l, 1)). Recall that in constructing standard MR models 
with RWA we have taken the spin operators for N a number of atoms as generators of su(2) , while for 
the radiation field we have considered bosonic realizations of different algebraic structures in (5). For 
EMG (9) however apart from su(2), only other possible algebra is su(l, 1), as given in (12). And even 
this can not be used, since differnt algebras can not be mixed in integrable Gaudin model. On the 
other hand, if we choose su(2) algebra for all operators and use the Holstein-Primakoff transformation 
(HPT) = tfyJf-N, S = yjf - Nb, S$ = N-S£, N = b^b for bosonizing the radiative mode, 
though we can retain the integrability, but get unphysical Hamiltonian for < \N\ > > 4^, with S^ 1 
becoming nonhermitian! 

4.2 High spin limit 

Fortunately, the above difficulties can be bypassed again by considering further a high spin limit 
so — » oo for the radiation field, which reduces the HPT to 

s ° + = i 6, ' s °" = ^ s « = -i?' e = ^°< (13) 

by retaining terms up to order O(^). 

For deriving now the integrable matter-radiation models without RWA , we replace spin operators 
at n = by bosons for the radiative mode in (11) by using (13) and choose the arbitrary parameters 
at n = 0,j as 

z = K + ex , zj = exj, j = 1, . . . N a , (14) 
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where K is the elliptic integral [23] related to a period of the elliptic functions. At the limit e — > 
therefore the coupling parameters w a (z n — z m ) can be expanded as 

k' 

w±(z - z k = -w±(z k - zq) -> ±— + O(e), w 3 (£o - z k ) = -w 3 (z k - zq) -»■ 1 + 0(e 2 ), 
- j( — ~ ~ ) + °( e )> w -( z j ~ z k) 0(e), u>s(«j - z fc ) -( — ~ - ) + O(e), (15) 

Denoting w± = \{w\ iti^), which reduce the elliptic Gaudin Hamiltonian to H n = -^Hn^ + \Hn^ + 
O(e ) with Hq^ = —J2jHj a \a = 1,2. Therefore, i^j a ^ give mutually commuting independent 
conserved quantities: [f/j a ' ) ,if^] = 0, with 



H?=sl H^=H^ + Hf S \j = l,...,N a , (16) 

where 

H? s) = n((bsl + tfs.) + ^sl + bsr)), n = ±=, (17) 

with elliptic modulus k = y/l — k' 2 , describes matter-radiation interaction having explicit CRW 
terms, where by an allowed phase transformation e~ in both b, SJ, we have changed the sign of the 
RW term. In (16) the part 



N a 1 

Hf s) = E ^^~M s k + s 7 s t + (is) 



accouts for interatomic interactions with inhomogenious coupling constants expressed through Xj. 
These arbitrary parameters may be adjusted to make the interaction strengths diminishing with 
increasing distance between the atoms. Various combinations of Hamiltonians (16) with arbitrary 
coupling constants: Hj = ujjSj,+Jj(Hj bS ^ + Hj SSS) ) can yield different integrable multi-atom matter- 
radiation models without RWA. For example, by taking simply J2jHj with loj = u,Jj = 1, we get 
an integrable multi-atom JC-type model with nontrivial CRW terms: 

#o = E + n ( b + 6t )(^ + + s l)) ( 19 ) 

3 

which in contrast to the popular belief [11] is exactly solvable through Bethe ansatz, as we see 
below. To make the model more physical we should add a field term Uftfb to the Hamiltonian (19). 
Considering the field frequency ojj to be small we can treat this additional term perturbatively over 
the Bethe solvable original integrable model. It is curious to note that, essentially at this small 
field frequency limit the RWA becomes a bad approximation and the CRW terms, which we have 
considered, become significant. 

5 Exact Bethe ansatz solutions 

In the basic algebraic Bethe ansatz (ABA) method the trace of the monodromy matrix r(A) = trT(X) 
produces conserved operators, while the off-diagonal elements ?2i(A) = B(X) and Ti2(A) = C(A) act 
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like creation and annihilation operators of pseudoparticles, with the M-particle state given as 
\M > B = B(Xi) ■ ■ ■ B(\m)\0 >, with the pseudovacuum |0 > defined through C(A)|0 >= 0. The 
major aim of ABA [19] is to find the eigenvalue solution: t(X)\M >b= A(A, {X a })\M 

All our models with RWA, constructed in a unified way through the ancestor model scheme, can 
be solved in a unified manner following the ABA method in its standard formulation [19]. Omitting 
here the details [17], we focus on the ABA solution of our integrable models without RWA (16- 19). 
For deriving the ABA relations of these models we have to repeat the steps we have followed in 
obtaining their Hamiltonians. That is we have to start from the ABA relations for the xyz spin 
chain with arbitrary spins , where a gauge transformed version of operators and states has to be used 
to define a modified pseudovacuum [23, 24]. Considering carefully the limit a — > the corresponding 
ABA relations for the elliptic Gaudin model is obtained [22]. Taking the high spin limit for the 
radiative mode in these relations we derive finally the eigenvalues and the Bethe equations for our 
integrable models without RWA. 

We start therefore from the ellipic Gaudin model, and using the key ABA relations from [22], we 
derive the energy spectrum as 

with arbitrary spin values s n of operators at the n-th site and involving elliptic function On. The 
corresponding Bethe equations can be found in [22] as 

E fln(Aq - z n ) u (X a - Xb) . , 

_ s "MA„- 2 „) + g e n (K - A S) (21) 

Choosing the inhomogeneity parameters at n = 0,j as in (14) , introducing a scaling also for the 
Bethe parameters: A a = el a and using some properties of the 9\\ function, we can derive now from 
(20) for n = j the exact eigenvalue for our matter-radiation model (16) without RWA, at the high 
spin limit (or equivalently, at e — ► 0) for the radiation field as 

I M -j N a a " \ 



where we have used short-hand notations #ii(0) = ofl, for I = 0,1,2. Similarly expanding (21) we 
extract the corresponding Bethe equation 

a" N a 1 M , 

^- w+ £*t^w-SpF«- (23) 

the solutions of which should determine the Bethe momentums lf,,b = 1, . . . , M for a given arbitrary 
nondegenerate set of inhomogeneity parameters Xj,j = l,...,N a . Note that xq can be absorbed by 
shifting the parameters Xj and l a . 

Either expanding (20) for n = and combining terms of different orders, or by taking the sum 
of (22) with an additional term with co and using the Bethe equation (23), we can obtain the exact 
eigenvalues for our multi-atom model (19), with explicit CRW terms as 

, N a' a" ( M N a \ 

E ( crw) = ^ (J2(x - k) + E s k (x -x k ))+ Wui/ (24) 



'10 



\b=l k=l 
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6 Concluding remarks 



In real matter-radiation models involving many atoms interatomic interactions must appear. Like- 
wise the customary rotating wave approximation in general must be avoided, at least when away 
from the resonance point. However these important cases could not be included without spoiling the 
integrability in most of the solvable models proposed earlier. We have proposed new integrable mod- 
els, where both these cases can be incorporated retaining the integrability of the system and extract 
exact solutions through the Bethe ansatz. For models with rotating wave approximation we find 
a series of integrable multi-atom models with inter-atomic interactions, which gives generalizations 
of Jaynes-Cummings, Buck-Sukumar and trapped ion models together with their q-deformations. 
These models belong to the rational or the trigonometric class and can be solved exactly in a unified 
way. 

We also construct for the first time exactly solvable multi-atom matter-radiation models with- 
out rotating wave approximation (RWA), which are close to the physical systems with the required 
counter rotating wave (CRW) terms. We achieve this by taking the high spin limit for the radiation 
field in the elliptic Gaudin model. However in contrast to our models with RWA, where we could 
construct different types of matter-radiation models, our success with models without RWA is re- 
stricted, where we could obtain only JC type models with CRW terms. Moreover, the free radiation 
term does not naturally appear in such integrable models and has to be included perturbatively over 
the exact BA result. 

Identifying the models in real systems and experimental verification of the related results pre- 
sented here, especially in many-atom microlasers [25] or in trapped ions away from resonance point 
[11] would be an important application. 
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